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Abstract. Previous smaller-scale studies of the globular cluster system of NGC 4636, an elliptical galaxy in the southern part 
of the Virgo cluster, have revealed an unusually rich globular cluster system. We re-investigate the cluster system of NGC 4636 
with wide-held Washington photometry. The globular cluster luminosity function can be followed roughly 1 mag beyond the 
turn-over magnitude found atV = 23.31±0.13 for the blue cluster sub-population. This corresponds to a distance modulus of 
(m — M) = 31.24 ± 0.17, 0.4 mag larger than the distance determined from surface brightness fluctuations. The high specific 
frequency is confirmed, yet the exact value remains uncertain because of the uncertain distance: it varies between 5.6 ± 1.2 
and 8.9 ± 1.2. The globular cluster system has a clearly bimodal color distribution. The color peak positions show no radial 
dependence and are in good agreement with the values found for other galaxies studied in the same filter system. However, 
a luminosity dependence is found: brighter clusters with an “intermediate” color exist. The clusters exhibit a shallow radial 
distribution within 7', represented by a power-law with an exponent of -1.4. Within the same radial interval, the galaxy light 
has a distinctly steeper profile. Because of the difference in the cluster and light distribution the specific frequency increases 
considerably with radius. At 1' and 9' the density profiles of the red and blue clusters, respectively, change strongly: the power- 
law indices decrease to around -5 and become similar to the galaxy profile. This steep profile indicates that we reach the 
outer rim of the cluster system at approximately 1 V. This interpretation is supported by the fact that in particular the density 
distribution of the blue cluster population can be well fit by the projection of a truncated power-law model with a core. This 
feature is seen for the first time in a globular cluster system. While the radial distribution of the cluster and field populations are 
rather different, this is not true for the ellipticity of the system: the elongation as well as the position angle of the cluster system 
agree well with the galaxy light. We compare the radial distribution of globular clusters with the light profiles for a sample 
of elliptical galaxies. The difference observed in NGC 4636 is typical of an elliptical galaxy of this luminosity. The intrinsic 
specific frequency of the blue population is considerably larger than that of the red one. * 

Key words, galaxies: elliptical and lenticular,cD - galaxies: individual: NGC 4636 - galaxies: star clusters - galaxies: stellar 
content - galaxies: structure 


1. Introduction 

The question as to what extent the study of globular cluster 
systems (GCSs) can contribute to the understanding of galaxy 
formation has been studied by many authors (e.g. Beasley et 
al. 120021 Cote et al. 120021 Forbes et al. 119971 Ashman & Zepf 
119921 but a satisfactory answer can still not be given. A guide 
to the answer could be the relation between global properties 
of GCSs and properties of their host galaxy. However, even 
when considering seemingly well-established statements, one 
faces complications when new and better observations become 
available. 

An example is the specific frequency Sn = N ■ 
1 0°-4(Mv—i5) introduced by Harris & van den Bergh (1 981 1 . 
with N being the total number of globular clusters (GCs) and 
My the absolute visual magnitude of the host galaxy. Theory 
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states that central galaxies in galaxy clusters should exhibit Sn- 
values about a factor 3 higher than galaxies in the field, which 
typically have values of 3-4 (e.g. Harris 120011 Elmegreen 1999 * 
Ashman & Zepf 1998 ). This mainly refers to 3 central galaxies, 
NGC 4486 (M87) in Virgo, NGC 1399 in Fornax and NGC 
3311 in Hydra I. But recent wide-held photometry of the GCS 
of NGC 1399 revealed only Sn ~ 6 (Dirsch et al. 2003a) , 
because the galaxy’s total brightness was previously underes¬ 
timated (confirming an earlier claim by Ostrov et al. 119981 . 
Although still somewhat higher than the average for ellipticals, 
it turns out to be significantly lower than in previous studies. 

NGC 4636, the southernmost bright elliptical galaxy in the 
Virgo cluster, is a counter example to the rule that high spe¬ 
cific frequencies are related to central cluster galaxies: it is 
a relatively isolated field galaxy and shows a peculiarly high 
Sn- In that and in other aspects, this elliptical galaxy is in¬ 
triguing. Because of its X-ray brightness, NGC 4636 has been 
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studied extensively in this frequency range and these observa¬ 
tions have also been used to study the dark matter distribu¬ 
tion. A surprisingly massive dark halo has been deduced by 
Loewenstein & Mushotzky ( 120031 : even well inside the effec¬ 
tive radius the dark matter constitutes a large if not the major 
fraction of the total mass, which is a unique case. Its bright, 
diffuse X-ray halo is atypical for an isolated galaxy and com¬ 
parable in luminosity and extension to a central cluster galaxy 
(e.g. Xu et al. 120021 Jones et al. 120021 Matsushita et al. 119981 
Awaki et al. [T994l . 

NGC 4636 also shows an unusual behavior in the far in¬ 
frared. Temi et al. 420031 ) found that its FIR luminosity be¬ 
tween 40 /rm and 80 //m is a factor of 50 higher than what 
is expected from the dust emission produced by its stellar 
body, a model that fits other elliptical galaxies quite well. Temi 
et al. speculate that a recently merged, dusty, gas-rich dwarf 
galaxy might explain the excess luminosity. Unfortunately, ex¬ 
isting HI observations cannot be used to support or rule out 
this idea: while Knapp et al. £[978}, Gallagher fl978l and 
Bottinelli & Gougenheim C1 9781 reported an extended HI disk 
with a mass of « 10 9 Mq, Krishna Kumar & Thonnard ( 1983 1) 
found only an upper limit of 10 8 Mq in NGC 4636. 

The globular cluster system of NGC 4636 has been inves¬ 
tigated by Hanes et al. fii)77l and Kissler et al. (H994t . They 
found very large specific frequencies of 9.9 and 7.5 ± 2.5, re¬ 
spectively. These values are typical for central galaxies and not 
for an elliptical in a low density region such as southern border 
of the Virgo cluster. With such a specific frequency NGC 4636 
would have one of the highest Sn values known. However, the 
study by Hanes et al. was based on photographic material and 
Kissler et al. only used V-band observations. It is therefore of 
great interest to confirm or revise the anomalous specific fre¬ 
quency by new deep photometry in more than one color band. 
A further objective is to determine the distance using the glob¬ 
ular cluster luminosity function (GCLF) and the spatial distri¬ 
bution of the GCs. 

As illustrated for NGC 1399 by Ostrov et al. (1 998 1 and 
Dirsch et al. (2 003 at. the specific frequency can easily be over¬ 
estimated due to an underestimation of the galaxy luminosity. 
Additionally, the specific frequency is known to to vary radially 
(e.g. Ashman & Zcpf 1998 Larsen et al. 120011 Rhode &Zepf 
120011 1. Hence it is necessary to observe the whole radial extent 
of a galaxy to derive this quantity. To reliably assess the specific 
frequency, it is thus crucial to obtain accurate measurements to 
large radii for both, the number of GCs and the light profile. 

Wide-held CCD studies exist only for five elliptical galax¬ 
ies: for two GC-rich central ellipticals: NGC 4472 (Rhode 
& Zenf l200H and NGC 1399 (Dirsch et al. 12003at and for 
three more “typical” ellipticals: NGC 3379, NGC 4406 and 
NGC 4594 (Rhode & Zenf 120031 . The studies of the central 
cluster galaxies revealed extended cluster systems that reach 
out to approximately 100 kpc. While the clusters in the inner 
part (less than roughly 50 kpc) have a distinctly shallower dis¬ 
tribution than the held population, there are indications in both 
cases that held stars and clusters are similarly distributed at 
larger radii. 

In this work we investigate the GCS of NGC 4636 for 
the first time with wide-held CCD images that cover roughly 


0.5° x 0.5° using the metallicity-sensitive Washington hlter 
system. This investigation aims to improve the understand¬ 
ing of the structural parameters of the GCS and its population 
structure. 

2. The data 

2.1. Observations & reduction 

The data set consists of Washington wide-held images obtained 
with the MOSAIC camera mounted at the prime focus of the 
CTIO 4-m Blanco telescope during 4 and 5 April 2002. We 
obtained four 600 sec images in R and seven 900 sec images in 
C. 

We used the Kron-Cousins R and Washington C filters, 
although the genuine Washington system uses T1 instead of 
R. However, Geisler (1996) has shown that the Kron-Cousins 
R hlter is much more efficient than T1 and that R and 
T1 magnitudes are closely related, with only a very small 
color term and zero-point difference. The MOSAIC wide-held 
camera images a held of 36' x 36'. The eight SITe CCDs 
have a pixel scale of 0.27"/pixel. Further information on the 
MOSAIC camera can be found on the MOSAIC homepage 
( http://www.noao.edu/kpno/mosaic/mosaic.htmh. 

We dithered the images to hll the gaps between the (eight) 
individual CCD chips. Due to the dithering not all of the area is 
covered by the same number of exposures, which restricts the 
usable held to 34(7 x 34(7. The seeing on the hnal images is 
1 ". 

The MOSAIC data has been handled using the mscred 
package within IRAF. In particular this software corrects for 
the variable pixel scale across the CCD which would otherwise 
cause a 4% variability of the brightness of star-like objects from 
the center to the corners. The hathelding resulted in images that 
had remaining sensitivity variations < 2.5% (peak-to-peak). 

To facilitate the search for point sources the extended 
galaxy light was subtracted. This was done using a median hl¬ 
ter with an inner radius of 9.5" and an outer radius of 11". This 
size is large enough to not alter the point source photometry 
which has been verihed with artificial star tests described later 
on. 

2.2. Photometry 

The photometry has been done using the PSF htting routine 
allstar within DAOPhot II. 

For the point source selection we used the \ and sharpness 
values from the PSF ht and found approximately 10000 point 
sources. The brightest non-saturated objects have T1 « 18.5 
(depending slightly on the individual MOSAIC chip). 

In each night 4 to 5 fields, each containing about 10 stan¬ 
dard stars from the list of Geisler (19961), were observed with a 
large coverage of airmasses (typically from 1.0 to 1.9). It was 
possible to use a single transformation for both nights, since the 
coefficients were indistinguishable within the uncertainties. 
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Fig. 1. The average completeness function of the entire held. 

We derived the following relations between instrumental 
and standard magnitudes: 

T1 =r + (0.71 ± 0.01) - (0.07 ± 0.01)X 
+ (0.033 ±0.003)(C-T1) 

(C - Tl) =(c - r) - (0.74 ± 0.02) - (0.20 ± 0.01)A 
+ (0.088 ±0.004)(C-T1) 

The standard deviation of the difference between instru¬ 
mental and calibrated magnitudes is 0.021 in Tl and 0.023 in 
C-Tl. 

A table containing the cluster candidates with coordinates 
and photometric data is available in electronic form. 

2.3. Photometric completeness 

The completeness of the data has been studied with the aid of 
the task addstar within DAOPhot II, which was used to add 
6666 stars to the science image. This number was found not 
to increase the crowding in the held. This was done ten times 
to produce ten different images. These modihed images were 
reduced in the same way as the original data. The hnal, spa¬ 
tially averaged completeness function for the entire MOSAIC 
held in the color range 0.9 < (C — Tl) < 2.5 is plotted in 
Fig.Qfor the entire MOSAIC held. The difference between the 
completeness function for red and blue clusters is marginal and 
does not need to be considered. However, there are strong spa¬ 
tial variations: close to the center of the elliptical galaxy and 
in the vicinity of bright foreground stars the completeness is 
lower. Throughout our study we applied completeness correc¬ 
tions as a function of spatial range. 

2.4. The color-magnitude diagram 

The color magnitude diagrams for all point sources in the 
MOSAIC images are plotted in Fig.|3 The left panel shows all 
objects within a galactocentric distance of 13'5 which we con¬ 
sider to be the appropriate distance to separate clusters from 
background objects (see Sect. 5.1). The so defined background 
objects are plotted in the right panel. The globular cluster pop¬ 
ulation is discernible within 0.9 < (C — Tl) < 2.1. The bi- 
modal nature of the color distribution (see Sect.3) is already 
visible. Background galaxies become signihcant at R> 23 and 
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Fig. 2. Left panel: the color magnitude diagram (CMD) of ob¬ 
jects with a radial distance smaller than 13(5. Right panel: the 
CMD of those objects further away than 13(5 from NGC 4636. 
We argue in Sect. 3 that in the more distant sample the number 
of globular clusters is negligible. This sample is used as back¬ 
ground. The globular cluster population is discernible within 
0.9 < (C-Tl) < 2.1. 

begin to dominate the CMD at R> 24. Most foreground stars 
are redder than the color range used for these CMDs. 

For the foreground absorption Schlegel et al. ( 119981 give 
Ab = 0.118 and Burstein & Heiles ( 1982 ) Ab = 0.05. Using 
the reddening law of Rieke & Lebofsky ( 1985 ) and Ec-ti = 
1.97 E b _v (Harris & Canterna 119771 . we find for the extinc¬ 
tion and absorption Ati = 0.08/0.04, Eq-ti = 0.06/0.03 
(Schlegel et al./Burstein& Heiles values), respectively. In the 
following we adopt the mean values, Ati = 0.06 ± 0.03 and 
Ec —ti = 0.04 ± 0.02. 

3. Cluster color distribution 

The color distribution of the cluster candidates is clearly bi- 
modal as shown in Fig.0 In this figure, as in the whole 
paper, the data is shown as histograms and as an adap¬ 
tive kernel filtered distribution (using an Epanechikov kernel, 
Merritt & Tremblay I1994F Whenever a fit is performed, the 
histogram data is used. This figure also shows how the dis¬ 
tribution varies with radius: the fraction of red clusters (those 
with 1.55 < (C — Tl) < 2.1) is higher at smaller radii 
(3(6 < r < 8(1) than at larger (8(1 < r < 13(5). A sum 
of two Gaussians is fitted to the color histograms: 

a ie -(( c - T1 )- b i) 2 /(2^) + a2e -((C-Tl)-b 2 ) 2 /(2<r 2 2 ). 

The fit results are tabulated in Tabled 

This fit shows that the peak position and the width of the 
Gaussians for the two populations remain constant with radius, 
despite the change of the relative heights. The radial constancy 
of the peak position has also been observed in the V-I color for 
a sample of elliptical galaxies by Larsen et al. (1200TV 
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radial range 



blue peak 




red peak 


3(6 < r < 8(1 

ai = 84 ± 4 

bi 

= 1.28 ±0.02 

in =0.15 ±0.02 

a 2 = 67 ± 4 

b 2 

= 1.77 ±0.02 

cr 2 = 0.19 ±0.02 

8(1 < r < 13(5 

ai = 46 ± 6 

bi 

= 1.28 ±0.02 

cri = 0.14 ±0.02 

a 2 = 16 ± 4 

b 2 

= 1.78 ±0.10 

cr 2 =0.20 ±0.10 


Table 1. The results of a two-Gaussian fit to the color distribution of the GC candidates brighter than R=24. 




C-T1 


C-T1 


Fig. 3. The color distribution of cluster candidates (brighter 
than R=24) is shown for an inner (3(6 < r < 8.' 1) and an 
outer (8(1 < r < 13(5) sample with solid and dashed lines, 
respectively. The sum of two Gaussians is fit to both distribu¬ 
tions, plotted with dotted lines. 

In the subsequent analysis we separate the clusters into two 
sub-samples, a red one (1.55 < (C — Tl) < 2.1) and a blue 
one (0.9 < (C — Tl) < 1.55) divided at the minimum of the 
inner sample plotted in Fig.0 

The color distribution depends on the luminosity: the frac¬ 
tion of bright red clusters compared to blue ones is larger than 
the corresponding fraction for faint clusters, indicating that the 
luminosity functions are different. This is shown in Fig0 The 
fitted peak positions are 1.28 ± 0.02/1.84 ± 0.04 for the faint 
sample and 1.31 ±0.04/1.69 ±0.02 for the bright sample. The 
color peak of the bright red clusters is approximately 0.1 mag 
bluer than that of the fainter cluster sample. We regard this as 
an indication that the concept of only two populations might 
not be sufficient. 

4. The luminosity function of the cluster 
population 

The LF of the cluster candidates is shown in Fig.0 for two 
radial intervals (2(3 < r < 5.'3 and 5.'3 < r < 11(2) and 
for the red (0.9 < (C — Tl) < 1.55) and the blue (1.55 < 
(C — Tl) < 2.1) cluster samples. The GCLF can be described 
by a Gaussian or a t5 function which has more extended wings 
(Seeker & Harris 179931 . However, no systematic differences in 
the turn over magnitudes (TOMs) have been reported when us¬ 
ing these functions (Larsen et al. 120011 1. We adopt the Gaussian 
model. The results of fitting the blue, the red and the entire 
cluster sample are shown in the first three columns of Tabled 


Fig. 4. The color distribution of cluster candidates shown for a 
bright (20 < Tl < 22) and a faint (22 < Tl < 24) sample 
with solid and dashed line, respectively. For both samples the 
sum of two Gaussians is fit, shown by the dotted lines. 


Since the position of the starting bin is arbitrary, we fitted two 
histograms with starting bins that are half the bin size apart and 
averaged the results which always deviated by less than 0.2 
mag. 

We have shown in Section 3 that the brighter clusters have 
a bluer red peak than the fainter ones. This indicates that two 
populations are not sufficient to describe the cluster system. 
We therefore divide the red cluster sample further, into an ’’in¬ 
termediate” color sample (1.55 < (C — Tl) < 1.8) and a 
“very red” sample (1.8 < (C — Tl) < 2.2). The ’’interme¬ 
diate” color sample selects the red peak of the bright cluster 
population. The results of a Gaussian fit to the data are found 
in Table[2] 

While most sub-samples are found to have similar TOMs 
and widths (which will be detailed shortly) the fit parame¬ 
ters for the intermediate color sample deviate significantly: the 
TOM is about one magnitude fainter and the width a factor 
of two broader. Moreover, the uncertainties are much larger. 
Despite the faint TOM we do not claim that there are relatively 
more faint clusters. We rather consider the faint TOMs with 
the enormous errors (at least a factor of 3 larger than in the 
other cases) as an indication that the fit of a Gaussian is not 
adequate. This interpretation is illustrated with Fig.[6]that com¬ 
pares the LF of the intermediate and the very red sub-samples. 
This figure shows that the fraction of bright GCs is larger for 
the intermediate color sample, consistent with the results using 
the color distributions. 

It is interesting to see this behavior in the context of 
the luminosity dependent color distribution of the GCS in 
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Fig. 6. The luminosity distributions of cluster candidates at 
larger radii for a red (1.8 < (C — Tl) < 2.15) and an ’’in¬ 
termediate” color (1.55 < (C — Tl) < 1.8) sample are plotted 
with solid and dashed lines, respectively. 


Fig. 5. The luminosity functions for the blue (0.9 < 

(C - Tl) < 1.55) and the red (1.55 < (C - Tl) < 2.1) GCs 
for an inner sample (2.'3 < r < 5.'3) in the upper panel and for 
an outer sample (5(3 < r < 11(2) in the lower panel. For all 
samples we plotted a histogram and a kernel smoothed distribu¬ 
tion. The errors of the histogram counts can be found in Table 2 
and were omitted in the figure. The average 50% completeness 
limit is marked by the vertical dotted line. 


NGC 1399: in this galaxy the bimodality disappears for the 
brightest GCs, which means a larger fraction of “intermedi¬ 
ate” colored GC candidates. The color peak of the brightest 
clusters in NGC 1399 is at C — Tl = 1.6 ± 0.1 while the 
“intermediate” colored sample of NGC 4636 is in the range 
1.55 < C — Tl < 1.8, hence the majority is redder than in 
NGC 1399. A further discussion of the potential nature of these 
sources is given in Sect. 7.2. 

Apart from the intermediate color sample no systematic dif¬ 
ferences in the TOM values can be seen within the uncertain¬ 
ties. This is in contrast to the expected dependence of the TOM 
on the metallicity of the cluster sample predicted by Ashman 
et al. ( 1995 ). This discrepancy is discussed in more detail in 
Sect. 7.2. The TOM that will be used from now on is the mean 
value of the blue (0.9 < (C — Tl) < 1.55) and the very 
red (1.8 < (C — Tl) < 2.1) samples in the entire field: 
Tl = 23.31 ± 0.13 and the width ct = 1.03 ± 0.15. 

The TOM and the width of the fitted Gaussian are strongly 
correlated and we find that TOMxi = (23.31 ± 0.13) + (0.7 ± 
0.11) • Act is a good fit where Act is the difference from the 
adopted width of 1.03. For example, for a width of ct = 1.3 
which is frequently kept constant when fitting the TOM of an 
elliptical galaxy in R which is comparable to Tl, the TOM 
would be Tl = 23.49 ±0.13. 


4.1. The distance to NGC 4636 

The absolute value of the TOM is best established in the V sys¬ 
tem (e.g. Richtler 2003] Harris l200Tl Ashman & Zenf 179981 . 
However, the V-R color of our clusters is unknown and the 
only possibility is to shift their R-TOM by a constant value. 
We assume a mean V — Tl = 0.50 ± 0.05 for our clus¬ 
ters (which has been derived from Galactic clusters of similar 
C — Tl color from the lists in Harris & Canterna ( !19771 > and 
Harris fT996t . Thus we find for the reddening corrected TOM 
V = 23.76 ±0.13 

Theoretically the TOM should depend on the metallicity 
(Ashman et al. 179921 However, in our sample no such ef¬ 
fect is visible. For this reason we do not attempt to correct 
for any metallicity dependence. According to Richtler ( 2003 1 
the TOM of the Milky Way, M31 and the mean of a (selected) 
sample of elliptical galaxies with SBF measurements is the 
same within the uncertainties (of the order 0.2 mag for each 
method), despite the possible metallicity differences between 
the GCS of these galaxies. He gives an averaged absolute TOM 
of —7.48 ± 0.11 in V. We hence find a distance modulus of 
(m-M) =31.24 ±0.17. 

Alternatively, we determined the R-TOM for the Milky 
Way using the GC compilation of Harris ((1 996 ) and find 
TOMr = —8.0 ± 0.16. With this value we determine a dis¬ 
tance modulus of (m—M) = 31.31±0.21, which is compatible 
with the above quoted value. 

The obtained distance modulus agrees with that of Kissler 
et al. ( 19941 . (m — M) = 31.4 ± 0.3 which is also based on 
the GCLF. However, the comparison with the distance modu¬ 
lus based on the SBF method is less favorable ((m — M) = 
30.83 ± 0.13, Tonry et al. l200H . The latter would correspond 
to a TOMti = 22.87 ± 0.18. We comment on this difference 
in the discussion. 
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radial 

interval 

(C — Tl) 

[0.9-2.1] 

(C — Tl) 
[0.9-1.55] 

(C — Tl) 
[1.55-2.1] 

(C - Tl) 

[1.55-1.8] 

(C - Tl) 

[1.8-2.15] 

1(3 - 11.'3 

23.38 ±0.14 

23.37 ±0.11 

23.41 ±0.21 

24.15 ±0.77 

23.24 ±0.14 


1.22 ±0.13 

1.17 ±0.12 

1.20 ±0.20 

2.10 ±0.50 

0.86 ±0.17 

iO 

1 

CO 

1—1 

23.44 ±0.18 

23.32 ±0.11 

23.63 ±0.27 

24.70 ± 1.30 

23.41 ±0.21 


1.26 ±0.16 

1.11 ±0.12 

1.44 ±0.23 

2.40 ± 0.90 

1.01 ±0.21 

CO 

1 

iO> 

23.36 ± 0.20 

23.42 ±0.25 

23.29 ±0.25 

24.02 ± 1.02 

23.18 ±0.15 


1.14 ±0.21 

1.14 ±0.23 

1.14 ±0.27 

1.94 ± 0.78 

0.79 ±0.17 


Table 2. The LF of GC candidates with T1 < 24 has been fitted with a Gaussian of variable width and TOM for several 
sub-samples. For each radial interval two lines are given, the first showing the TOM, the second the width of the Gaussian (a). 
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Fig. 7. The radial distribution of all cluster candidates, the blue 
(0.9 < (C - Tl) < 1.55) and the red (1.55 < (C - Tl) < 
2.1) sub-samples are shown in the three panels from top to bot¬ 
tom, respectively. The limiting magnitude is Tl=24. We also 
included the radial cluster density profile of Kissler et al. ( 1994 1 
in the uppermost panel using cross-symbols. Their data has a 
similar limiting magnitude and the data points can be compared 
directly. The galaxy light profile in Tl (in arbitrarily scaled in¬ 
tensities, see Sect. 5.2) is shown as dashed line in the upper 
panel. Power-law fits are shown as straight dashed lines. Open 
circles mark the point source density profile uncorrected for 
background. The backgrounds are indicated by the horizontal 
dotted lines. 

5. Spatial distribution 

5.1. Radial cluster distribution 

To determine the radial density distribution we accounted for 
the radially variable completeness by constructing the LF of 
all point sources brighter than Tl=24 and correcting the LF 
for completeness. The LF is then summed, normalized to the 
area and background corrected. Beyond r > 13(5 the number 
density of the point sources in the whole/blue/red color inter- 


r [arcmin] 

all clusters 

blue clusters 

red clusters 

1' -8' 

8' - 12' 

1' - 3(5 
3(5-6' 

6' - 12' 

-1.44 ±0.08 

-4.9 ±0.3 

-1.52 ±0.04 
-5.1 ±0.3 

-0.78 ±0.04 
-1.73 ±0.14 
-4.9 ±0.4 


Table 3. This table lists the exponents of power-law fits 
(n(r) ex r") to the cluster surface density determined for dif¬ 
ferent radial intervals. The blue clusters can be described by 
two exponents, while three exponents are needed for the red 
cluster and the complete sample. 


val remains constant and hence it has been chosen as back¬ 
ground. The used background densities are 2.33 ± 0.07/1.88 ± 
0.06/0.76 ± 0.04[obj/arcmin 2 ] for the entire/blue/red sample. 
The resulting radial profiles are shown in Fig. 0 and are tabu¬ 
lated in Table lA.3l In addition, the values uncorrected for back¬ 
ground and the background values are plotted in Fig. |7] which 
illustrates the uncertainties of the employed backgrounds. 

In all radial distributions shown in Fig.Qa change in slope 
at a radial distance between 7' and 9' is apparent. This steep¬ 
ening occurs approximately two arc-minutes closer in for the 
red than for the blue sub-population. It is already visible for 
the blue population in the uncorrected sample, however not as 
clearly. As it is customary we fit power-laws to the data which 
are shown in Fig.[7]as dashed lines. To account for the variable 
exponents that manifest themselves in the different slopes in 
the double logarithmic presentation we split the radial range in 
smaller parts. The resulting exponents (n(r) cx r Q ) are tab¬ 
ulated in Tabled While the blue clusters can be described 
with one exponent within 1' — 8' the density distribution of 
the red ones changes the exponent within this radial interval. 
The uncertainties of the power-law indices are given for the 
background used (including the statistical error). However, a 
systematic under- or overestimation of the background density 
might significantly change the slopes for the larger radius inter¬ 
val. Since an underestimation is not very probable considering 
Fig-El the main uncertainty is an overestimation, i.e. the cluster 
system is more extended. However, beyond 13' a radial depen¬ 
dence is not visible in any of the sub-samples nor does the color 
distribution of the point sources show the characteristic peaks 
of the GCS. 

In the uppermost panel of Fig.[7]we also included the den¬ 
sity profile by Kissler et al. ( 1994 1. They used V measurements 
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with a limiting magnitude of 24. Since V — R « 0.5 (Galactic 
globular clusters with C-T1=1.5) this corresponds roughly to 
our brightness limit. Except for the last point the agreement is 
excellent. The last point deviates since they apparently overes¬ 
timated their background because they used a too small CCD 
field for a reliable background determination. We also included 
in Fig. □ the galaxy light profile measured in the T1 filter (see 
the next section). 

The observed density profile can be described with the pro¬ 
jection of a truncated power-law model with a core and a uni¬ 
form exponent a (truncation radius r t , core radius r c ). Within 
the core a constant number density is assumed. The resulting 
distributions are shown in Fig. □for the red and the blue pop¬ 
ulations. In particular for the blue population such a simple 
model provides an excellent fit to the observations. The (used) 
parameters for the blue/red populations are: a = —2.3/ — 2, 
r c = l.'l/H5 and r t = 12710'. 

Finally, one might ask whether the ’’intermediate” red sam¬ 
ple (1.55 < (C — Tl) < 1.8) has the same distribution as the 
red and blue sample. Its radial distribution (a = —1.5 ± 0.2, 
r < 7') is indistinguishable within the uncertainties from these 
other distributions. 


5.2. Galaxy profile 

The surroundings of NGC 4636 are devoid of any large galaxy 
or bright foreground stars and it is straightforward to use the 
ellipse task in IRAF/STSDAS to model the galaxy. In Fig.0 
the luminosity and color profiles are shown and tabulated in 
Table lA.4l Caon et al. 1 1994) presented a B luminosity pro¬ 
file for this galaxy which we include in our figure. The V-Tl 
color resulting from the surface profile of Fauer et al. <19951 is 
V-T1= 0.60 ± 0.05, which is in excellent agreement with the 
color given in Prugniel et al.’s ( 199 8) compilation (0.58±0.02). 
This agreement supports the photometric calibration of Tl. The 
same holds within r < Of5 for the B-Tl color that results from 
Caon et al.’s profile. However, for larger radii Caon et al.’s pro¬ 
file is slightly shallower and the B-Tl color becomes bluer than 
C-Tl. However, this behavior is not supported by the color pro¬ 
file of Idiart et al. ( 20021 . 

A good empirical fit to the Tl profile is 


Tl(r) 


—2.5 log 





with ai = 3.3 • 10 -7 , a 2 = 5.5 • 10" 9 , ri = 0.11', r 2 = 8.5', 
cti = 2.2, o 2 = 7.5. The standard deviation between the fit 
and the Tl profile is 0.037 and the maximal deviations less 
than 0.09 mag (r < 8'). 

To further check our photometric calibration we compare 
the integrated galaxy light with aperture measurements com¬ 
piled by Prugniel et al. <19981 . shown in Fig.^ For this com¬ 
parison we transformed our Tl measurements into R. The over¬ 
all agreement is very good. However, for radii larger than 0.'5, 
the measurements from the literature seem to be 0.03 magni¬ 
tudes fainter than ours. Since this shift is within our calibration 
uncertainty we make no attempt to correct for it. 



0.01 0.05 0.50 5.00 

major axis [arcmin] 



major axis [arcmin] 

Fig. 8. In the upper panel we show the galaxy luminosity profile 
measured in Tl (solid line) and C (dashed line). The dotted line 
shows the B profile determined by Caon et al. (11994 ) and the 
dash-dotted line the V profile from Fauer et al. (Q995). In the 
lower panel the resulting color profiles are shown. A solid line 
is used for C-Tl and a dotted one for B-Tl. 



Fig. 9. The integrated luminosity profile is compared to aper¬ 
ture measurements collected by Prugniel et al. ( 1998 ). 


In Fig. Hit is obvious that for radii smaller than 7' the clus¬ 
ter population has a distinctly shallower distribution than the 
diffuse light, a fact that has already been noted for the in¬ 
ner 4' by Kissler et al. ( 1994 ). At larger radii, the slope of 
the galaxy light and cluster density profile appear very simi¬ 
lar; however, the uncertainties in the background determination 
make a firm statement difficult. Beyond 11' (which corresponds 
to « 60 kpc) the stellar population of NGC 4636 is negligible 
especially compared to the X-ray halo which reaches 300 kpc 
(Matsushita et al. fT998l . 
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The color gradient of NGC 4636 is relatively shallow (for 
a comparison with a galaxy that has a strong radial color gra¬ 
dient see for example NGC 1399, Dirsch et al. 12003al l. Such 
a gradient can be caused by differences in the distribution of 
distinct stellar populations, however, in NGC 4636 high reso¬ 
lution archival WFPC2 HST images reveal an abundance of 
patchy dust structures in the inner region. Therefore, reddening 
can play a major role in producing the color gradient. Hence, 
the stellar populations seem to be well mixed. 

5.3. Azimuthal distribution 

To determine the ellipticity of the cluster system we study 
the azimuthal density distribution of GC candidates in cir¬ 
cular annuli. The ellipticity e and the number density along 
the major and the minor axis (N ffl ,Nh) are related via e = 
1 — (Nb/Na) 1 ^ 0 , where a is the exponent of the radial den¬ 
sity distribution (r~ a ). 

In Fig.[j0l the azimuthal cluster density is shown for four 
different radial bins. In such a graph an elliptical cluster dis¬ 
tribution manifests itself as a sin(20 + (f>o) distribution. We fit 
this function to the data which is shown in Fit’. llOl as the dot¬ 
ted line. Within 2(3 and 7'9 the distribution is clearly elliptical. 
Within a radial range of 7(9 and 10(7 a fit also results in an 
elliptical distribution, however, the distribution is considerably 
asymmetric and the fit is relatively poor. At even larger radii the 
distribution appears rather asymmetric with an enhancement of 
clusters towards the North-West (around 310°). This excess can 
be interpreted as a slightly more extended cluster distribution 
towards this direction. Interestingly, the galaxy light also shows 
an indication for such an asymmetric elongation, however, flat- 
field uncertainties prohibit a definite answer, see Fig. El 

The radial dependence of the ellipticity and the position an¬ 
gle is shown with filled circles in Fig^] Within 8' the cluster 
system has an ellipticity comparable to that of the galaxy light 
(e ss 0.35). At larger radii no information for the galaxy light is 
available and the cluster ellipticity drops to 0.23 ± 0.09, only a 
la deviation from the ellipticity found at for the smaller radii. 

It is also of interest to study whether a difference in the az¬ 
imuthal distribution of red and blue clusters can be found. For 
this task we used all blue/red clusters within 2(3 and 8', the 
radial range where the two samples are indistinguishable. We 
found that the values for ellipticity and position angle agree 
with the values obtained for the total sample. Also the ’’inter¬ 
mediate” color sample defined above has the same azimuthal 
distribution within the uncertainties. 

6. The total number of clusters and the specific 
frequency 

To determine the total number of clusters we used the TOM 
and the width of the Gaussian that had been fit to the GC lu¬ 
minosity function in Sect. 4 (TOM(Tl) = 23.31). With this 
model we extrapolate from the number of observed GCs to 
all clusters. We proceeded as follows: we summed all clusters 
that are brighter than Tl=24 in radial bins. These cluster counts 
were completeness corrected to account for the radially chang¬ 
ing completeness. In this way we observed 66% of the total 



0 50 100 150 200 250 300 

4 > 


Fig. 10. The azimuthal distribution of the cluster candidates 
(0.9 < (C - Tl) < 2.1, 20 < T1 < 24) is plotted for four dif¬ 
ferent radial bins (2(3-5(1,5(1-7(9,7(9-10(7,10(7-13(5). 
The dotted line shows the fit of a function asin(20 + b) + c, 
which describes the number counts of an elliptically distributed 
population that are counted in circular bins. The position angle 
is counted from North to East. 


number of clusters. The final number of all clusters within a 
radial distance of 14' is 4200 ± 120. 

Kissler et al. ( 11994l i found 3600 ± 500 within 7', where we 
find 3500 ± 170, which is in excellent agreement. 

The specific frequency is the number of globular clusters 
per unit luminosity scaled to a galaxy with an absolute lumi¬ 
nosity ofMy = — 15mag(SN = N-10°- 4 ^ Mv+15 ^). In Sect.5.2 
we have shown that V—Tl = 0.60±0.05. The light profile can 
be observed to a semi-major axis distance of 9(5. The GCs are 
counted in circular rings and for the comparison we need the 
galaxy light in circular rings as well. The largest aperture for 
which we can determine a reliable luminosity is 7(8. At this ra¬ 
dius the total Tl luminosity of NGC 4636 is Tl = 8.70 ±0.05. 
This is the maximum radius for which a specific frequency can 
be obtained. The total number of cluster candidates within this 
radius is 3570 ± 180. Hence we find Sn = 5.8 ± 1.2 and 
Sn = 8.9 ±1.2 for the GCS and the SBF distance, respec¬ 
tively. A slight inconsistency in these values is that we used the 
same number of clusters and varied only the galaxy’s luminos¬ 
ity. However, since the TOM is visible the number of clusters 
does not change for the shorter distance. 

Since the clusters have a shallower distribution than the dif¬ 
fuse galaxy light the specific frequency varies radially and it 
is not sufficient to give only one number for the specific fre¬ 
quency. Therefore we calculated the specific frequency for dif- 
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r [armin] 

Fig. 11. In this figure we compare in the upper panel the el- 
lipticity of the cluster system, measured in three radial bins, 
with the ellipticity of the galaxy light determined in the R im¬ 
age (solid line), the C image (dashed line) and in a B image by 
Caon et al. C 19941 dotted line). In the lower panel the same is 
done for the position angle. As x-axis we used the geometrical 
mean of semi-major and semi-minor axis. 


r[arcmin] 

Sn 

(m- M) = 31.25 

Sn 

(m - M) = 30.83 

0.5 

3.7 ±0.2 

5.7 ±0.3 

1.4 

3.5 ±0.2 

5.4 ±0.3 

2.3 

3.9 ±0.2 

5.9 ±0.3 

3.2 

4.3 ±0.2 

6.5 ±0.3 

4.1 

4.6 ±0.2 

7.1 ±0.3 

5.0 

4.9 ±0.2 

7.5 ±0.3 

5.9 

5.3 ±0.3 

8.0 ±0.4 

6.8 

5.6 ±0.3 

8.5 ±0.4 

7.7 

5.8 ±0.3 

8.9 ±0.4 


Table 4. Specific frequency based on a V-R=0.6 color without 
the the distance uncertainty that is constant ASn = 1.2. 


ferent limiting radii for the distance of Tonry et al. and 

ours. The result is shown in Fin. II H and tabulated in Tabled 
The specific frequency is rather large, in particular if one 
uses the smaller distance of Tonry et al. ( 2001 1. In this case 
NGC 4636 would have one of the largest specific frequencies 
ever observed. In our opinion this is a further argument against 
this shorter distance. 


Fig. 12. The R image of NGC 4636 is shown (33' x 33'), North 
is up. East to the left. The galaxy light appears to be more ex¬ 
tended towards the North-West. 



Fig. 13. The specific frequency is plotted for different limiting 
radii. Open circles are used for a distance modulus of 30.83 
(Tonry et al. 120011 and filled circles for (m — M) = 31.25 
(GCLF). The errors are only statistical errors. The distance un¬ 
certainties result in an additional uncertainty of ASn = 1-2. 

context. While definitive answers cannot be given, we attempt 
to work out the environment inside which possible solutions 
can be searched for. We summarize our results in Tabled 

7. 1. The specific frequency 


7. Discussion 

The properties of the GCS of NGC 4636 present several puz¬ 
zles. In the following discussion we put them into a broader 


One of these puzzles is the specific frequency, which is surpris¬ 
ingly high for a relatively isolated galaxy, even if our large dis¬ 
tance modulus of 31.25 is correct. A distance modulus of 30.8, 
as suggested by the SBF method, shifts Sn to unrivaled val- 
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Coordinates (J2000): 

12“ 

42 m 49.9 s , 2°41'16" 

Assumed reddening: 

Ec- 

_xi = 0.04 ±0.02 

TOM distance modulus: 

(m M ) = 31.24 ±0.17 

SBF distance modulus: 

(m M ) = 30.83 ± 0.13 

Total number of GCs: 

N - 

= 4200 ± 120 

Total apparent magnitude: 

T1 

= 8.70 ±0.05 

Average color: 

V- 

- T1 = 0.60 ±0.05 

Global specific frequency: 



TOM distance: 

Sn 

= 5.8 ± 1.2 

SBF distance: 

Sn 

= 8.9 ± 1.2 


Table 5. Basic data for NGC 4636 used from the literature or 
derived in this work 


ues. Richtler ( 200311 compared SBF distances and GCLF dis¬ 
tances for a large sample of early-type galaxies. Apart from 
an overall good agreement, deviations are almost exclusively 
found in the sense of smaller SBF distances. These galaxies 
frequently exhibit other indicators pointing to the existence of 
intermediate-age or even young stellar populations. The differ¬ 
ence of 0.4 mag between SBF and GCLF moduli, as found for 
NGC 4636, is only marginally explainable by the quoted mea¬ 
surement uncertainties. Moreover, the appearance of a super¬ 
nova (Giclas 119391 is another indication for an intermediate- 
age population. On the other hand, in those galaxies where 
younger populations have unambiguously been identified, for 
instance in NGC 1316 (Goudfrooij et al. l200~fl or NGC 5018 
(Hilker & Kissler-Patig 119961 Leonardi & Worthey [2000 1 , the 
specific frequency is distinctly lower due to the higher lumi¬ 
nosity of younger populations (Gomez & Richtler [20011 ). It is 
therefore difficult to explain the difference between SBF and 
GCLF moduli by a contribution of an intermediate-age popu¬ 
lation to the SBF magnitude. A small contribution would not 
be sufficient given the insensitivity of the I-band fluctuation 
magnitude to population differences (e.g Cantiello et aI.[2003}, 
while a large one is improbable because of the high specific 
frequency. The specific frequency of NGC 4636 is comparable 
to that of NGC 1399, the central galaxy in the Fornax cluster 
(Dirsch et al. l2003al . This demonstrates that a central position 
in a high density environment is not a necessary condition for 
developing a high ,Sly-value. Scenarios like stripping globu¬ 
lar clusters from neighboring galaxies, which have been con¬ 
structed to explain the wealth of clusters for example in NGC 
1399, are not applicable in the case of NGC 4636. 

The extent and brightness of the X-ray halo as well as the 
total mass indicated by the X-ray analysis could also be char¬ 
acteristic of a small group of galaxies. Together with the high 
specific frequency, one may ask whether NGC 4636 is the rem¬ 
nant of a galaxy group. Isolated elliptical galaxies have sev¬ 
eral times been suspected to result from multiple merging of 
small groups (e.g. Stocke et al. 120041 Mulchaey & Zabludoff 
119991 Vikhlinin et al. 119991 . Numerical simulations of these 
merger processes show a strong dependence on the initial con¬ 
ditions (Athanassoula et al. 119971 , suggesting that the merger 
rate for group galaxies encompassed by a common dark matter 
halo is low, the corresponding timescales possibly exceeding a 
Hubble time. Thus NGC 4636 might be the merger outcome 
of a close group of galaxies with individual dark halos and the 


present halo the dispersed debris of the merger components. 
More properties, which support such a scenario, are the unusu¬ 
ally high dust content (Temi et al. 120031 and the absence of 
a color gradient which suggests that chemical evolution was 
largely terminated at the formation time of NGC 4636. 

The high, perhaps extremely high, star formation rate that 
accompanied the merger is favorable for a high efficiency 
of cluster formation (Larsen & Richtler l2000l . giving a qual¬ 
itative explanation for the high specific frequency of NGC 
4636. Searching for possible counterparts, one comes across 
NGC 1132, which has been suspected of being the merger rem¬ 
nant of a small galaxy group (Mulchaey & Zabludoff 119991 . 
However, nothing is known about its cluster system due to its 
large distance of about 100 Mpc. 

McLaughlin (1 9991 ) argued that the efficiency of GC forma¬ 
tion is constant if one normalizes the total cluster mass to the 
sum of stellar mass and X-ray gas mass. Does the high specific 
frequency come back to ’’normal” when the gas mass contained 
in the X-ray halo of NGC 4636 is included? For this, one needs 
factors of the order 1.5-2. One problem is the uncertain volume, 
inside which this normalization has to take place. The GCS ter¬ 
minates at a radius of roughly 40 kpc. According to Matsushita 
et al. ( 1998 k the stellar mass inside this radius still is a factor 
of 10 larger than the X-ray gas mass, rendering the inclusion 
of the gas mass insignificant. Beyond this radius, the galaxy’s 
light cannot be traced any more. Only at a radius of 200 kpc, 
apparently much larger than the extension of the GCS and the 
galaxy, the X-ray gas mass becomes comparable to the total 
stellar mass. The question remains how the history of GC for¬ 
mation and the X-ray gas are related. A more or less constant 
ratio of GC mass and total mass does not contradict a local high 
efficiency of GC formation in a merger event. 

7.2. The luminosity function 

It has been theoretically predicted by Ashman et al. ( I1995fe that, 
for a given break in the mass function, metal-poor clusters have 
a brighter TOM than more metal-rich clusters. Such differences 
have also been observed in the V luminosity function of sev¬ 
eral GCS by Larsen et al. ( 2001 ). Despite those expectations 
we do not see convincing evidence for this effect in NGC 4636. 
For example, considering the blue and the very red population 
(see Table[2} and using the Washington color-metallicity rela¬ 
tion of Harris & Harris ( 200 21. we find a mean metallicity of 
-1.44dex and -0.08dex for these populations. Such a metal¬ 
licity difference should manifest itself in a TOM difference 
of ATI « 0.3 (we used the fact that T1 is basically R and 
Table 3 of Ashman et al. fT995ii . Such a difference cannot be 
seen in Tabled However, the fitted width of the Gaussian is 
much smaller for the very red than for the blue sample. If we 
had used the same width we would have found a TOM for the 
whole radial range of T1 = 23.46 ± 0.18 which is also within 
the uncertainties equal to that of the blue sample. We have no 
good explanation for the absence of the predicted TOM depen¬ 
dence on metallicity. 

We have shown that the intermediate color cluster sample 
has a LF that differs from that of the other sub-populations 
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in the sense that it has brighter GCs (Fig. [3. This behavior 
can also be seen in the brightness selected color distributions 
(Fig.0. It is interesting to note that this is similar to the GCS 
around NGC 1399, where the brighter GCs have a different 
color distribution (Dirsch et al. 12003al Ostrov et al. 11998 ;). 
In particular the bright GCs have no bimodal color distribu¬ 
tion any longer but peak at an intermediate color. We tenta¬ 
tively interpret both observations with a picture in which the 
“normal” GCs are mixed with a population of objects that are 
brighter and of intermediate color. These may be former nuclei 
of stripped nucleated dwarf galaxies. 

7.3. The’’edge" of NGC4636 

The GC density profile becomes very steep beyond 8' (see 
Fig-0- Since the GCS traces the galaxy light (see also the 
Section7.4), it is reasonable to assume that the whole stellar 
body of NGC 4636 shows this behavior, which is supported by 
a relatively sharp border of the galaxy light. Such steepening in 
the very outer regions of elliptical galaxies has been explained 
by Jaffe (1 987 ) and White ( 11987 1 on the basis of energy con¬ 
siderations. They argue that the transition from an r -3 to an 
r -4 spatial density profile occurs if there is a sharp cut in the 
energy distribution function near the escape energy in the case 
of an approximately Keplerian potential. 

A different case is, for example, NGC 1399, which has also 
been studied with wide-held CCD photometry by Dirsch et al. 
( 2003a ). They attributed its (spatially) uniform r -3 profile to 
its extended dark matter halo. The Keplerian regime (if it exists 
at all in the inner region of the Fornax cluster) would then be 
beyond the observable radial range. If this idea was applicable 
to NGC 4636 it would imply that we have already reached the 
regime of a declining potential. 

However, this interpretation conflicts with the existence of 
a massive dark halo - as is claimed on the basis of X-ray ob¬ 
servations (Loewenstein & Mushotzkv 120031 Matsushita et al. 
119981 Mushotzky et al. 119941 . Hence either our interpretation 
of the light profile and of the cluster distribution is incorrect 
or the X-ray mass profile is inaccurate. Problems concerning 
the interpretation of the X-ray gas arise from its turbulent na¬ 
ture (Otho et al. 120031 Jones et al. 1 2002 1 and a possible inter¬ 
action with the Virgo inter-cluster gas (Trinchieri et al. 119941 ). 
Therefore an independent measurement of the dark matter dis¬ 
tribution is desirable. 

7.4. Is the GC radial density distribution peculiar? 

The radial distribution of the GCS in NGC 4636 is consid¬ 
erably shallower than the galaxy light over nearly the entire 
radial range. How can this be assessed in a quantitative man¬ 
ner? GCSs generally exhibit shallower surface density profiles 
than their host galaxies’ light (Harris 120011 Ashman &Zepf 
fT998i> . It is also known that larger galaxies tend to have 
GCSs with shallower radial distributions (Kissler-Patig JT997| 
Ashman & Zepf 199 81. which also holds for the galaxy profile 
(Kormendy & Diorgovski tl 989t . However, we are not aware of 
any recent quantitative comparison examining these trends. In 


the following we collect power-law indices fitted to the radial 
distribution of GCSs from the literature and perform a simi¬ 
lar, homogeneous fit to published galaxy light profiles within 
a comparable radial range. The correlation of these power-law 
indices with galaxy luminosity is determined and compared for 
GCSs and field stars. 

We use the galaxy profile measurements of Virgo and 
Fornax ellipticals published by Caon et al. (1 994 ) and fit power- 
laws to the data (we multiplied the luminosity profiles by -0.4 
to ensure comparability with the cluster data). A power-law 
fit may not be the preferred fit to a galaxy light profile. It is 
chosen here to allow a comparison with the published clus¬ 
ter density profiles which are generally described by power- 
laws. In order to exclude galaxies for which a single power- 
law at radii larger than 20" is not appropriate, we fitted two 
radial ranges and omitted those galaxies for which the power- 
law indices deviate by more than 0.15. The B luminosity pro¬ 
files from Caon et al. d 19941 were transformed to V magni¬ 
tudes with B-V colors from Prugniel et al. ( 1998 1. We adopted 
the total apparent magnitudes given by Caon et al.. For the 
absolute magnitudes we used the distance moduli of Tonry 
et al. ( 2001 ). For galaxies without a SBF distance we used a 
mean cluster distance modulus of 31.4 for Fornax and 31.1 for 
Virgo. To correct for absorption the mean of the Schlegel et 
al. (i 19981 ) and Burstein & Heiles ( 1982 ) values was employed. 
The result is shown in the upper panel of Fig. ED which il¬ 
lustrates the well-known fact that larger early-type galaxies 
are less concentrated (e.g. Kormendy & Diorgovski ll989l ). In 
this plot we also included for a consistency check the same 
dependence published by Schombert et al. 119861) . This data, 
however, is much more heterogeneous and not reddening cor¬ 
rected. The distances are those provided by Schombert et al. 
( 1986 ) and are based on the recession velocities. Furthermore, 
the galaxy luminosities published by Schombert et al. are mea¬ 
sured within a linear radius of 16kpc. To transform these lu¬ 
minosities to the scale of Caon et al. we compared galaxies 
that are contained in both samples and obtained a mean shift 
of my Caon = my Schombert — 0.45. One expects a systematic 
trend of this shift with magnitude, however, we were unable to 
detect it, presumably due to the small magnitude coverage of 
the galaxies in common (12 objects). We fit a linear function to 
the galaxies brighter than My = —20 of the two samples and 
found a slope of —0.31 ± 0.03, consistent with Schombert’s 
sample. Considering its luminosity (My = —21.78 ± 0.13), 
NGC 4636 has a surprisingly shallow light distribution (a = 
-1.63 ±0.05). 

To compare the dependence of the GCS radial distribution 
and the luminosity profile on the total galaxy luminosity we 
compiled data from the literature, the most important single 
source being Kissler-Patig ffW7l) in conjunction with distance 
moduli from Tonry et al. ( 2001 1. Furthermore, we excluded 
those cases in which the error of the power-law exponent is 
larger than 0.3 and cases in which the power-law fit is done on 
the basis of less than 4 points and included a few galaxies with 
recent measurements. For the galaxy luminosity we have cho¬ 
sen the value given by Caon et al. ( 1994 1 or, when not available, 
the total V magnitude from the RC3. The resulting dependence 
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of the cluster radial distribution on the host galaxy’s luminosity 
is shown in the lower panel of Fig. M 

The comparison of the radial cluster distribution with the 
galaxy light shows the well established fact that in general the 
cluster systems have shallower distributions. However, it also 
shows that the difference between the slopes of galaxy light and 
the GC radial profiles depends on luminosity: brighter galaxies 
exhibit a slightly larger difference on the average, which was 
not yet obvious in the sample discussed by Ashman & Zepf 
FT9981 

The radial distribution of the GCs of one elliptical is rather 
different: NGC 1316, the brightest galaxy in the sample, clearly 
stands out. However, it is not a typical elliptical galaxy, but a 
merger remnant with an age of 2-3 Gyr (Goudfrooij et al. l200ll 
Gomez & Richtler 12001> . 

The location of NGC 4636 in this diagram on the other hand 
is typical. 

In Tablc lA.iil we also included the difference between the 
slope of the GC radial density profile and the galaxy light for 
several galaxies where both numbers are available. It shows 
that both slopes are closely related, irrespective of luminosity. 

It is important to note that a shallower density profile of the 
GCS does not necessarily imply a GCS that is more extended 
than the galaxy light. Both in NGC 4636 and NGC 1399 the 
GCS has a clearly shallower profile than the galaxy light until 
a certain radius. At large radii, however, the distribution be¬ 
comes steeper and is compatible with that of the galaxy light. 
The radius at which this happens is fairly large and hence the 
shallower inner decline does not correspond to a core. 

7.5. The discrepancy between the distribution of 
clusters and field stars 

Another fact worth pointing out is that the ellipticities of 
the GCS and the field population of NGC 4636 are the same 
(within the uncertainties) within 4' and 9' while the radial dis¬ 
tributions are rather different in this range. 

Also puzzling is the fact that no pronounced radial color 
gradient can be seen for the galaxy light. One might expect that 
any field component that is connected to the metal-poor, blue 
cluster population (and which should therefore follow the same 
density profile) would become radially more prominent and 
hence the galaxy light is expected to have a stronger color gra¬ 
dient than actually observed. The implication is that the popula¬ 
tion that has formed together with the clusters (in particular the 
metal-poor, blue ones) is negligible compared to the field popu¬ 
lation that formed with a very low cluster formation efficiency. 
This is in strong contrast to other galaxies where the clusters 
have a shallower distribution and in which strong color gradi¬ 
ents in the light profile have been found: NGC 1427 (Forte et 
al. l200Tll . NGC 4472 (Rhode & Zepf l200TI >. NGC 1399 (Dirsch 
et al. l2003at . 

In the following we perform a study similar to that pre¬ 
sented in Dirsch et al. ( 2003al l: we aim at constraining the so- 
called intrinsic specific frequency (ISF) which measures the 
specific frequencies of the red clusters with respect to a red 
population of the same mean color and of the metal-poor clus- 



Mv 



-20 -21 -22 -23 -24 


Mv 

Fig. 14. Upper panel: the dependence of the power-law index 
of the surface brightness profile on galaxy luminosity (multi¬ 
plied by -0.4) of Virgo/Fornax galaxies published by Caon et 
al. (1 994!) (see text for a detailed description). For the galaxies 
plotted as filled circles distance determinations via SBF exist 
(Tonry et al. 120011 while mean cluster distances have been as¬ 
sumed for the galaxies plotted with open triangles. The solid 
line is a fit to these galaxies brighter than My = —20. The dots 
are similar measurements obtained by Schombert et al. C l 986 1. 
transformed to the same magnitude scale (see text). The dashed 
line shows the fit to these data (also brighter than My = —20). 
Lower panel: the exponents of power-law fits to the GC den¬ 
sity distributions of elliptical galaxies is plotted (filled circles) 
together with the galaxies from the upper sample having SBF 
distances (open circles). The crosses show the data of Kissler- 
Patig’s compilation with errors in the exponent larger than 0.3 
or in which the distribution has been obtained with less than 4 
data points. It is obvious that the GCs in general have a shal¬ 
lower distribution than the galaxy light and that the correla¬ 
tion with the total galaxy luminosity is very similar between 
them. The data plotted with the filled symbols is tabulated in 
Table lA.5l and Tab I c l A. til The solid and dashed lines show the 
fits to the galaxy light and correspond to the lines in the upper 
panel. The dotted line shows a fit to the GCS radial distribution. 


ters with respect to their underlying field population. Such an 
approach results in a unique solution only when the two popu¬ 
lations are assumed and the two observables, light profile and 
color profile, are fitted by varying only the blue and red ISF. 
The result is shown in Fig. El For radii larger than 6' a con¬ 
stant ISF for each population can be used, which corresponds 
to the fact that cluster and light distribution become similar at 
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Fig. 15. Upper-most panel: the observed light profile (solid 
line) is compared to the fitted light profile (solid circles), based 
on an assumed ISF, which is different for the red and blue clus¬ 
ters (see text for details). Middle panel: the color profile re¬ 
sulting from the fit (solid circles) is shown together with the 
observed profile. Lower panel: the ISF of the blue and the red 
population used to generate the luminosity and color profile is 
shown. The ISF values are transformed into the V system using 
a constant V-R color of 0.6. This has been done for an easier 
comparison with the values published in Dirsch et al. ( 2003 a). 


large radii. The ISF of the blue population is, as expected, at all 
radii larger than that of the red ones. Furthermore, at large radii 
the value of the blue population’s ISF is very similar to that 
derived for NGC 1399. However, while in the latter galaxy the 
blue ISF remains constant with radius, it declines in NGC 4636 
inwards. This decline is accompanied by a decline in the red 
cluster density profile, otherwise a color gradient would arise. 
The value of the red ISF is considerably smaller in NGC 4636 
than in NGC 1399. 

Summarizing, we find that in particular the red clusters 
have formed with a lower efficiency in NGC 4636 compared 
to NGC 1399, while it is slightly higher for the blue clusters 
for radii greater than 6'. Within this radius also their formation 
efficiency was lower. Despite these differences the global spe¬ 
cific frequencies are very similar and hence the overall cluster 
formation efficiencies have been comparable. 

7.6. The globular cluster - X-ray gas connection 

The GCS and the hot X-ray gas may have a common evolu¬ 
tionary history (via mass loss during star formation and SN 
explosions of the respective stellar population), however, there 


is no reason to expect a tight relationship in their distribution. 
It is therefore striking to find that similar radii appear to be 
characteristic of both. We have shown than the GC density dis¬ 
tribution becomes steeper at 6' — 8' (depending on the popula¬ 
tion) which is also where the X-ray gas changes its slope and 
becomes flatter (Trinchieri et al. 119941 . At a radius of about 
13' the cluster system appears to have ended; this is the radius 
where the intensity profile of the X-ray gas changes again and 
becomes nearly flat (at ss 12kpc, Matsushita et al. 119981 . Is 
this merely a coincidence? A study of the dynamical properties 
of the GCS that can constrain the matter distribution is required 
to shed light on these findings. 

Appendix A: The data tables 

The luminosity function of the GCs, their color and radial dis¬ 
tribution and the galaxy profile including the ellipticity are pub¬ 
lished electronically. The galaxy sample used in Sec. 7.4 is also 
available electronically. The table captions are given in this ap¬ 
pendix. 

Table A.l. The position, the T1 magnitude and C-Tl color for 
each of the point sources identified in the investigated frame. 
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C-Tl 

N(3(6 < r < 8(1) 

N(8.'l < r < 13(5) 

0.89 

4.2 ±5.8 

-7.9 ±9.7 

0.96 

9.7 ±6.5 

10.5 ± 11.0 

1.02 

7.0 ± 7.4 

-3.1 ±9.7 

1.10 

23.4 ±7.7 

29.2 ± 12.4 

1.17 

63.8 ±10.0 

31.1 ± 12.4 

1.24 

83.6 ±11.0 

35.4 ± 12.8 

1.31 

75.8 ± 10.5 

49.1 ± 13.1 

1.38 

79.0 ± 10.5 

46.3 ± 12.3 

1.45 

61.7 ±9.4 

17.5 ± 10.4 

1.52 

39.6 ± 7.7 

16.5 ±9.3 

1.58 

40.8 ± 7.9 

9.7 ±9.2 

1.66 

61.5 ±9.0 

33.4 ± 10.0 

1.73 

73.2 ±9.7 

-0.4 ±8.5 

1.80 

63.7 ±8.7 

15.1 ±7.6 

1.87 

56.2 ±8.4 

13.6 ±8.1 

1.94 

32.6 ±6.4 

16.5 ±6.8 

2.00 

36.9 ±6.9 

9.8 ±6.8 

2.08 

11.3 ±4.6 

6.8 ±6.5 

2.15 

7.3 ±4.2 

-1.2 ±5.8 

2.21 

10.0 ±4.3 

-0.6 ±5.3 


Table A.2. The color distribution of the GC candidates brighter 
than R=24 is shown in Fig.0 In this table we used the his¬ 
togram counts with a bin-width of 0.07 mag. 
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r 

n(all)[arcmin z ] 

with background background subtracted 

n(blue)[arcmin 

with background background subtracted 

n(red)[arcmin 

with background background subtracted 

0'45 

64.25 ± 3.39 

61.93 ±3.37 

21.63 ±2.15 

20.87 ±2.15 

50.04 ± 2.92 

48.16 ±2.92 

1'35 

38.59 ± 1.47 

36.26 ± 1.47 

14.43 ± 0.87 

13.67 ±0.88 

26.81 ± 1.25 

24.94 ± 1.26 

2'25 

23.20 ± 0.69 

20.87 ±0.69 

9.71 ±0.47 

8.95 ± 0.48 

14.77 ±0.54 

12.90 ± 0.55 

3'15 

16.70 ±0.48 

14.36 ± 0.49 

7.91 ± 0.33 

7.16 ±0.34 

9.40 ± 0.37 

7.52 ± 0.38 

4' 05 

12.26 ± 0.54 

9.93 ±0.55 

5.52 ± 0.35 

4.76 ± 0.35 

7.56 ± 0.47 

5.69 ±0.47 

4'95 

8.90 ± 0.40 

6.57 ±0.40 

3.56 ± 0.28 

2.80 ±0.28 

5.44 ± 0.29 

3.57 ±0.29 

5'85 

8.01 ±0.41 

5.69 ±0.41 

3.39 ±0.26 

2.63 ±0.26 

5.07 ± 0.34 

3.20 ±0.35 

6'75 

6.72 ±0.33 

4.39 ±0.34 

2.85 ±0.22 

2.09 ±0.22 

4.34 ± 0.28 

2.46 ±0.28 

7'65 

5.01 ±0.25 

2.68 ±0.26 

1.72 ±0.15 

0.97 ±0.15 

3.68 ±0.22 

1.81 ±0.23 

8'55 

4.68 ±0.23 

2.35 ±0.24 

1.45 ±0.13 

0.69 ±0.14 

3.54 ±0.20 

1.66 ±0.21 

9'45 

3.84 ±0.20 

1.51 ±0.21 

1.18 ±0.12 

0.42 ±0.12 

2.94 ±0.18 

1.07 ±0.19 

10'35 

3.30 ±0.21 

0.98 ±0.22 

1.10 ±0.12 

0.35 ±0.13 

2.48 ±0.18 

0.60 ±0.19 

11'25 

3.00 ±0.20 

0.65 ±0.19 

0.99 ± 0.09 

0.24 ±0.10 

2.23 ±0.16 

0.36 ±0.17 

12'15 

2.65 ±0.17 

0.32 ±0.18 

0.75 ±0.09 

0.00 ±0.10 

2.22 ±0.15 

0.34 ±0.16 

13'05 

2.85 ±0.18 

0.52 ±0.20 

0.93 ±0.11 

0.17 ±0.12 

2.35 ±0.16 

0.47 ±0.18 

13'95 

2.05 ±0.15 

-0.27 ±0.16 

0.70 ±0.10 

-0.05 ±0.11 

1.62 ±0.13 

-0.26 ±0.14 

14'85 

2.33 ±0.16 

0.01 ±0.17 

0.75 ±0.09 

0.00 ±0.10 

1.88 ±0.15 

0.00 ±0.16 

15'75 

3.15 ±0.25 

0.82 ±0.26 

0.85 ±0.14 

0.09 ±0.15 

2.77 ±0.22 

0.89 ±0.23 

16'65 

2.79 ± 0.23 

0.46 ±0.24 

0.94 ±0.13 

0.18 ±0.13 

2.05 ±0.20 

0.17 ±0.21 

17'55 

3.19 ±0.26 

0.86 ±0.27 

1.07 ±0.17 

0.32 ±0.17 

2.47 ±0.21 

0.60 ±0.22 

18'45 

2.78 ±0.31 

0.45 ±0.32 

1.18 ±0.22 

0.42 ±0.23 

1.87 ±0.25 

0.00 ±0.26 

19'35 

2.02 ±0.27 

-0.31 ±0.28 

0.36 ±0.14 

-0.40 ±0.14 

2.12 ±0.30 

0.24 ± 0.30 


Table A.3. The radial densities of all GCs and the red and blue subsamples with and without background subtraction as shown in 
Fig m The densities given here are calculated for clusters brighter Tl=24. To include the entire luminosity function these values 
have to be multiplied by 1.52 ( see Sect. 5). 
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r 

T1 

[mag/arcsec -2 ] 

C-Tl 

PA 

e 

(TOO 

16.27 ±0.03 

1.91 ±0.05 

88 ±9 

0.31 ±0.08 

0(01 

16.34 ±0.03 

1.91 ±0.05 

144 ± 12 

0.04 ±0.02 

0(02 

16.46 ±0.03 

1.88 ±0.05 

108 ±9 

0.03 ±0.01 

0(03 

16.76 ±0.03 

1.86 ±0.05 

124 ±5 

0.03 ±0.01 

0(05 

17.01 ±0.03 

1.85 ±0.05 

120 ± 20 

0.00 ±0.01 

0(07 

17.33 ±0.03 

1.81 ±0.05 

157 ±3 

0.02 ±0.01 

0(10 

17.63 ±0.03 

1.77 ±0.05 

155 ±3 

0.02 ±0.01 

0(14 

17.97 ±0.03 

1.75 ±0.05 

153 ±2 

0.03 ±0.01 

0(18 

18.35 ±0.03 

1.74 ±0.05 

153 ±2 

0.04 ±0.01 

0(20 

18.61 ±0.03 

1.74 ±0.05 

150 ±2 

0.04 ±0.01 

0(26 

18.89 ±0.03 

1.73 ±0.05 

151 ±2 

0.06 ±0.01 

0(34 

19.27 ±0.03 

1.69 ± 0.05 

150 ± 2 

0.08 ±0.01 

0(45 

19.65 ±0.03 

1.67 ± 0.05 

149 ±2 

0.11 ±0.01 

0(54 

19.94 ±0.03 

1.67 ± 0.05 

147 ±2 

0.13 ±0.01 

0(65 

20.23 ±0.03 

1.67 ±0.05 

145 ±2 

0.14 ±0.01 

0(77 

20.52 ±0.03 

1.66 ±0.05 

146 ±2 

0.17 ±0.01 

1(04 

20.89 ±0.03 

1.66 ±0.05 

145 ± 2 

0.22 ±0.01 

1(38 

21.28 ±0.03 

1.66 ± 0.05 

145 ±2 

0.25 ±0.01 

1(65 

21.61 ±0.03 

1.65 ± 0.05 

146 ±2 

0.25 ±0.01 

2(00 

21.94 ±0.03 

1.62 ± 0.05 

145 ±2 

0.27 ±0.01 

2(20 

22.19 ±0.03 

1.60 ±0.05 

145 ±2 

0.26 ±0.01 

2(56 

22.56 ±0.03 

1.59 ±0.05 

143 ±2 

0.26 ±0.01 

3(10 

22.92 ±0.03 

1.59 ±0.05 

144 ±2 

0.28 ±0.01 

3(41 

23.14 ±0.03 

1.59 ±0.06 

144 ±2 

0.28 ±0.01 

4(14 

23.50 ±0.03 

1.59 ±0.07 

146 ±2 

0.33 ±0.01 

4(54 

23.80 ±0.03 

1.53 ±0.10 

146 ±2 

0.33 ±0.01 

5(01 

24.07 ±0.03 

1.53 ±0.10 

146 ±2 

0.33 ±0.01 

5(49 

24.28 ± 0.03 

1.50 ±0.12 

145 ± 2 

0.35 ±0.01 

6(04 

24.63 ± 0.03 


144 ±2 

0.35 ±0.01 

6(66 

24.92 ± 0.03 


144 ±2 

0.35 ±0.01 

7(31 

25.19 ±0.04 


147 ±2 

0.35 ±0.01 

8(04 

25.41 ±0.05 


157 ±3 

0.37 ±0.02 

8(86 

25.49 ± 0.05 


156 ±4 

0.43 ±0.01 

9(73 

26.01 ±0.06 


156 ±4 

0.39 ±0.02 


Table A.4. The galaxy luminosity profile, color profile, posi¬ 
tion angle (PA) and ellipticity (e) is tabulated. The errors on 
the position angle and ellipticity are those given by the task 
ellipse. The luminosity and color uncertainties include calibra¬ 
tion errors as well. 


galaxy 

M v 

UGC 

Lit 

NGC 524 

-21.90 ±0.22 

-1.71 ±0.12 

KP97 

NGC 708 

-22.41 ±0.32 

-1.00 ±0.11 

OH02 

NGC 720 

-22.07 ±0.18 

-2.2 ±0.1 

KP96 

NGC 1052 

-21.09 ±0.30 

-2.26 ±0.27 

KP97 

NGC 1316 

-23.43 ±0.19 

-2.04 ±0.20 

GeaOl 

NGC 1374 

-20.53 ±0.15 

-1.80 ±0.20 

KP97 

NGC 1380 

-21.40 ±0.18 

-1.60 ±0.10 

KP97b 

NGC 1399 

-22.48 ±0.19 

-1.40 ±0.15 

Dea03a 

NGC 1404 

-21.67 ±0.23 

-2.00 ±0.10 

KP97 

NGC 1427 

-21.03 ±0.26 

-2.00 ±0.20 

FeaOl 

NGC 1700 

-22.54 ±0.25 

-1.10 ±0.07 

BeaOO 

NGC 3115 

-21.23 ±0.10 

-1.84 ±0.23 

KP97 

NGC 3258 

-21.45 ±0.30 

-1.60 ±0.20 

Dea03b 

NGC 3268 

-21.68 ±0.32 

-1.80 ±0.20 

Dea03b 

NGC 3379 

-20.93 ±0.11 

-1.73 ±0.15 

RZ03 1 

NGC 4278 

-21.01 ±0.24 

-1.85 ±0.20 

KP97 

NGC 4365 

-22.06 ±0.18 

-1.15 ±0.25 

KP97 

NGC 4406 

-22.41 ±0.14 

-1.24 ±0.05 

RZ03 

NGC 4472 

-23.21 ±0.10 

-1.28 ±0.04 

RZ01 

NGC 4486 

-22.54 ±0.16 

-1.61 ±0.08 

KP97 

NGC 4594 

-22.14 ±0.18 

-1.85 ±0.07 

RZ03 

NGC 4636 

-21.78 ±0.24 

-1.30 ±0.10 

this paper 

NGC 4697 

-21.30 ±0.18 

-1.90 ±0.20 

KP97 

NGC 7014 

-21.82 ±0.34 

-1.71 ±0.11 

OH02 

NGC 5193 

-21.99 

-1.38 ±0.07 

OH02 

IC 5193 

-23.49 

-1.02 ±0.08 

OH02 


Table A.5. Total absolute luminosities of the galaxies (taken 
from either the Caon et al. H994I - transformed into V - or 
from the RC3.9) and the power-law exponent of its GCS that 
are plotted as filled circles in the lower panel of Fig^] KP97 
- Kissler-Patig ( 1997 1. KP97b - Kissler-Patig et al. (i l997bt . 
KP96 - Kissler-Patig et al. Q996}, RZ01 - Rhode &Zepf 
(120011 RZ03 - Rhode &Zepf i i200l . Fea91 - Forte et al. 
(2001), Dea03a - Dirsch et al. ( 2003at . Dea03b - Dirsch et al. 
I 2“003b1 . GeaOl - Gomez et al. ( 2001 1 . BeaOO - Brown et al. 
d2000l l. OH02 - Okon& Harris (120021 . For the last two galax¬ 
ies no SBF distance is available and thus the distance used in 
the paper on the GCS was assumed. 1 We fitted the GC den¬ 
sity profile using the data points of RZ03 with radial distances 
smaller than 11'. 
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galaxy 

Mv 


dgal 


Ogal 

- 

QiGC 

NGC 1336 

-19.09 

± 

0.21 

-2.53 

± 

0.09 




NGC 1351 

-20.30 

± 

0.16 

-2.31 

± 

0.04 




NGC 1373 

-18.62 

± 

0.47 

-2.74 

± 

0.07 




NGC 1375 

-19.36 

± 

0.13 

-2.58 

± 

0.02 




NGC 1379 

-20.47 

± 

0.15 

-2.58 

± 

0.09 




NGC 1380 

-21.40 

± 

0.18 

-2.01 

± 

0.04 

0.39 

± 

0.11 

NGC 1380A 

-18.70 

± 

0.29 

-2.65 

± 

0.04 




NGC 1381 

-19.95 

± 

0.21 

-2.61 

± 

0.05 




NGC 1399 

-22.48 

± 

0.16 
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± 

0.01 

0.27 

± 

0.15 

NGC 1404 

-21.67 

± 

0.19 

-2.28 

± 

0.05 

0.28 

± 

0.10 

NGC 1419 
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± 

0.24 

-2.85 

± 

0.03 
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-21.03 

± 

0.24 

-2.22 

± 

0.04 

0.22 

± 

0.20 

NGC 4261 

-22.47 

± 

0.19 
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0.08 




NGC 4472 
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± 

0.06 

NGC 4600 

-17.04 

± 

0.22 

-2.40 

± 

0.07 




NGC 4636 

-21.78 

± 

0.13 

-1.62 

± 

0.05 

0.32 

± 

0.21 
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